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We establish relations between the statistics of g factors and the fluctuations of energy in metallic 
nanoparticles where spin-orbit coupling is present. These relations assume that the electron dynam- 
ics in the grain is chaotic. The expressions we provide connect the second moment of the g factor to 
the root-mean square "level velocity" (the derivative of the energy with respect to magnetic field) 
calculated at magnetic fields larger than a characteristic correlation field. Our predictions relate 
readily observable quantities and allow for a parameter-free comparison with experiments. 
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PACS numbers: 73.23.Hk, 71.70.Ej 

It was noted in experiments 1 " 2 that the Zeeman split- 
ting of discrete energy levels in nanoparticles is very 
sensitive to the presence of spin-orbit interaction. The 
splitting can be characterized by a level-dependent g fac- 
tor. Adding to the Al grains only 4% of Au resulted in 
a change of measured g factors from approximately 1.7 
to 0.7. In addition to being suppressed, the g factor in 
the presence of spin-orbit interaction also fluctuates ran- 
domly from level to level. 

These fluctuations of g factors were described in the 
framework of random matrix theoryS*^ (RMT) and the 
suppression of the g factor was related to the strength 
of the spin-orbit interaction and to the elastic mean free 
path of electrons in the grains— The limit of strong spin- 
orbit interaction corresponds to a short spin-orbit scat- 
tering time, T so 6/h <C 1, where S^ 1 is the mean den- 
sity of states in the grain at zero field. It was predicted 
that the distribution of eigenvalues of the g-factor ten- 
sor in this limit should have a Gaussian form. Within 
RMT, this distribution is characterized by a phenomeno- 
logical parameter (g 2 )- In Ref. |j, this parameter was 
expressed in terms of the grain size, electron mean free 
path I, and relaxation time r so . However, a comparison 
of the experimental results with theory was not entirely 
satisfactory. On the one hand, there is an indication that 
the distribution of the eigenvalues and directions of the 
eigenvectors of the g-tensor is Gaussian and correspond 
to a "pure" symplectic ensemble—^ On the other hand, 
it is not clear— whether the small values of (g 2 ) obtained 
in the experiments agree with the theoretical values esti- 
mated in Ref. 0. The difficulty in making the comparison 
comes from the lack of information about the amount of 
disorder in the grain. 

The goal of this work is to provide three relations for 
the distribution width (g 2 ^ to other quantities that are 
directly measured in the same set of experiments (and 
thus do not rely on any additional information about the 
amount of disorder in the grains). These quantities are 
the variance of the energy level derivative with respect 
to the magnetic field (known as level velocity^) and the 



zero-magnetic field level curvature (the second derivative 
of the energy level with respect to magnetic field). 

We begin by stating our main results. Our first expres- 
sion, valid for strong spin-orbit coupling only, r so 5/H <C 
1, is 
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where var[] denotes the variance, /is is the Bohr mag- 
neton, and B* is the crossover field for breaking time- 
reversal symmetry. Equation Q gives a statistical con- 
nection between the response at B — > (the g factor) 
to that at large magnetic fields (the level velocity). It 
provides a way to check experimentally if the grains ex- 
hibiting fluctuations of the g factor indeed belong to a 
"pure" symplectic ensemble, rather than to an ensem- 
ble describing the crossover between the orthogonal and 
symplectic limits. Provided that one collects data for the 
dispersion of energy levels over a sufficiently large range 
of magnetic fields, the terms on both sides of Eq. (JTJ can 
be found independently using the same data set. This 
expression is universal and contains no microscopic or 
materials parameters. 

The second expression we find relates properties of two 
sets of grains which are equivalent macroscopically except 
for the value of r so . It reads 
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where go denotes the materials bulk value for the g factor 
(.go = 2 for free electrons). In Eq. (s 2 ) is evaluated in 
the strong-spin-orbit coupling limit (r so — > 0). However, 
the second term on the right-hand side is evaluated for 
t so — > oo (absence of spin-orbit coupling). For example, 
one may consider two sets of Al:Au grains, one with no 
doping and another with moderate dopingi 

The two terms on the right-hand side of Eq. (J2J come 
from distinct contributions. The first term, the "spin 
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part," is associated with the Debye mechanism of energy 
dissipation associated with spin reorientations; the sec- 
ond one, the "orbital part," is due to the eddy currents 
induced in the grain. We note that the spin-orbit scat- 
tering rate 1 /r so apparently may be estimated for a given 
material or host-dopant pair. 6 Equation © permits us 
to separate the spin and orbital contributions to the fluc- 
tuations of the g factor. 

The second moment of the g factor can also be related 
to the statistics of other spectral quantities, such as the 
zero-field level curvature 
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As Eq. ©, this expression is applicable in the strong- 
spin-orbit coupling regime only (symplectic ensemble). 

It is important to remark that the g factor in a given 
metallic nanoparticle is in reality a tensor with three dif- 
ferent eigenvalues, even for grains that are statistically 
isotropic^*! However, the distribution of the matrix el- 
ements of the g-factor tensor is still characterized by a 
single quantity (g 2 )- Throughout our manuscript, when 
establishing relations between (g 2 ) to other quantities, 
we consider a magnetic field applied in some fixed but 
arbitrary direction. Once these relations are established, 
one may use Ref. to construct the full statistics of the 
g-factor tensor. 

We will now establish Eqs. ©-©. The main idea 
behind the derivations is to relate statistical quantities 
to invariants of the system, such as the traces of the 
magnetic moment operator. For that purpose, let us be- 
gin by writing the Hamiltonian for the disordered (or 
chaotic) grain in the presence of an applied magnetic field 
as H(B) = Ho + B M, where the magnetic moment oper- 
ator has both orbital and spin parts: M — M or b + Mspin- 
To simplify the discussion, we assume that the grain is 
isotropic. We define the g factor of the nth energy level 
as 
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where {e na } are the eigenvalues of Hq. Note that due 
to Kramers degeneracy at B = 0, the levels are twofold 
degenerate. We use the index a to distinguish states 
that are time-reversal conjugate. The matrix elements 
of Ho follow either the symplectic (/3 = 4) or orthogo- 
nal (/3 = 1) ensemble statistics, depending on whether 
spin-orbit coupling is present or absent, respectively. For 
both cases, the matrix elements of M, when expressed in 
the eigenbasis {|ncr;0)} of H , fluctuate according to a 
Gaussian distribution with zero mean. For the symplectic 
ensemble, the variance of the diagonal matrix elements 
reads 4 



\{na;0\M\na\Q)\' 



3Tr 
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where the trace runs over the 2N states in the band, 
with TV 3> 1 being assumed (the factor of 2 accounts 
for Kramers degeneracy). The degeneracy at zero field 
allows us to pick a basis such that (na; 0\M\na'; 0) is di- 
agonal in the a indices. Using Eq. © and first-order 
perturbation theory, we find that (na; 0\M\na'; 0) = 
(— 1)°" 8 aa i g n UB/2- Thus, from Eq. ©, we arrive al£ 
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Notice that the quantities on the left-hand side of Eqs. 
© and © are defined at B = 0. 

Now we have to write the statistical quantities that 
appear on the right-hand side of Eqs. ©-©, in terms of 
Tr(M 2 ). Note that the latter is an invariant and therefore 
takes the same value at zero or large magnetic fields. 

Let us first consider Eq. ©. The variance of the level 
velocity can be computed in terms of the variance of the 
matrix elements of the magnetic moment operator since 
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For a sufficiently large magnetic field Bo S> B*, time- 
reversal symmetry in the grain is broken. In the presence 
of strong-spin-orbit scattering, orbital and spin degrees of 
freedom remain mixed, but the ensemble statistics of the 
Hamiltonian eigenstates switches from 27V x 27V symplec- 
tic to 2N x 27V unitary ((3 = 2). Thus, we need to com- 
pute the variance of the matrix elements of the magnetic 
moment operator in the unitary regime. For this pur- 
pose, we make use of the eigenvalues and eigenvectors of 
the magnetic moment operator: M\ka) = (— l) a Mk \ka), 
with k = 1, . . . , TV and a = ±1 due to the time-reversal 
properties of This yields 



(v; B \M\v; Bo) = ^(-l) a Af fe | (v; B \ka) 
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For the unitary ensemble in the large- N limit, the eigen- 
vector amplitudes shown on the right-hand side of Eq. 
© fluctuate independently according to the Porter- 
Thomas distribution»ii One finds that |(^;Bo|^a)| 2 ) = 
and (|(^; Bo\ka)\ 4 ) = \N 2 , independently of state 

indices. Hence, the average matrix element of M must 
vanish and the variance can be written as 
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Putting together Eqs. ©, ©, and ©, we arrive at Eq. 
©. 

To derive Eq. © , we separate the magnetic moment in 
terms of spin and orbital parts, M = Af spin + Af or b , which 
are statistically independent from each other. From Eq. 
©, we obtain 
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The spin contribution can be written in terms of the 
imaginary part of the ac spin susceptibility of a free 
electron gas in the presence of spin-orbit coupling. 4 The 
susceptibility can then be evaluated using conventional 
means of rate equations at frequencies much larger than 
the mean level spacing in the grain, yet smaller than the 
spin-orbit scattering rate. In the limit of t so S/H <C 1 one 
finds^ 
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This corresponds to the first term on the right-hand side 
ofEq. ©. 

The orbital contribution in Eq. (|10|> is ensemble in- 
dependent and is the same regardless of the presence or 
absence of the spin-orbit interaction. It is convenient to 
evaluate it in the limit of r so —> oo to decouple spin and 

orbital degrees of freedom. To relate the Tr ^M„ T ^\ to 

the variance of the velocity of spin-resolved levels at large 
fields, we can use Eq. Q and note that M sp in just in- 
duces a constant slope in the dispersion of the energy 
levels with magnetic field at r so — > oo, while all fluctua- 
tions are caused by M OT h- Therefore, we write 
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B = Bn 



to isolate the fluctuating part (note that Bo 3> B* here 
as well). It is important to observe that in the absence of 
spin-orbit mixing and for large magnetic fields the statis- 
tics of the eigenstate |n Ti-Bo) corresponds to a N x N 
unitary ensemble (rather than to 2N x 2N when spin 
and orbital parts are strongly coupled) . Similarly to Eq. 
©, we can decompose the matrix elements in Eq. <|12[) 
using the eigenstates of M rb- In this case, however, 
due to the change in the ensemble dimension, we have 
(|(nT;Bo|fcA)| 2 ) = l/N and (| (n T; B \kX) | 4 ) - 2/N 2 
for large N, independently of the orbital quantum num- 
ber n. Finally, combining all these results, we arrive at 
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Inserting this expression into Eq. JUJ, we obtain the 
second term on the right-hand side of J2J. We remark 
that Eq. © is fully consistent with Eq. (31) of Ref. [H 
where a description of the statistical properties of the g 
factor including the intermediate crossover regime was 
developed in terms of phenomenological RMT parame- 
ters. 

In order to obtain Eq. we follow an approach simi- 
lar to that employed in the two previous derivations. We 
define the level curvature at B = as 
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(To simplify the notation, here we set e na = e n .) Since 
eigenvalues and eigenfunctions fluctuate independently in 
the Gaussian ensembles, we find that 
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with n n' and a, a' taking arbitrary values. The pref- 
actor in Eq. I|15fl is defined as 
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,(£n-e«') 2 , 

where the delta function is used to fix the energy level 
in the middle of the band. The average of eigenvalues 
can be performed using the appropriate two-level cluster 
function^ In the limit N — > oo, we find 

= w> (17) 

0=4, y 

with d denoting the inverse of the mean density of 
states^ The ensemble average of off-diagonal matrix el- 
ements of the magnetization is also easily computed in 
terms of the trace of the magnetization operator in the 
limit of large N: 
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and 
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Inserting Eqs. (|17|l -119 (l into (|15|l we arrive at 
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For the symplectic ensemble, one can showi^*i& that 
\J {K 2 ) = (n\/3/4:)(\K n \). Thus, using this relation and 
combining Eqs. © and we obtain Eq. 0. 

Equation J3}, like Eq. Q), also involves only quanti- 
ties that are directly measurable. However, in practice, 
the difficulty in obtaining large statistics for the second 
derivative at B = from the tunneling conductance data 
makes it less appealing when applied to experiments^! 

Once the variance of the level velocity is obtained 
from the experimental data, it may also allow for an- 
other test of RMT. Consider the level velocity correlation 
function 8 ^ 
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For a pure ensemble, this correlation function can be 
rescaled to a universal form. Defining the correlation 
field as 



B c = 1/y/CJfl) 



(22) 



and calling x = AB/B C and c(x) = B 2 C C V {AB), the di- 
mensionless correlation function in the unitary ensemble 
has the asymptotes^ 



c{x) 



1 - 2tt 2 x 2 , x -C 1, 
-l/(vrx) 2 , x > 1. 



(23) 



The full shape of the correlation function is known 
from numerical simulations* 8 *^ as well as from analyt- 
ical calculations^ 

Finally, it is interesting to note that the correlation 
field is related to the amount of disorder in the grains 
when the electron motion is diffusive. It is straightfor- 
ward to show that at r so — > oo, 



B r = K 



^>o/L 2 



(24) 



where $o is the flux quantum, L is the grain linear size, 
kp is the Fermi wavelength, and k is a dimensionless 
coefficient that depends on the grain geometry. For a 
spherical shape, k = 3n / (2^/5) , in which case L is the 
grain radius. By measuring the variance of the level ve- 
locity at large fields, one can obtain C„(0) and find the 



experimental value of B c from Eq. (|22|l . Using Eq. i|24[) . 
one can then get an independent estimate of the amount 
of disorder present in the grain. Another approach is to 
fit the universal curv o 8 i 18 c{x) to the experimental data 
and obtain C v (jS) as a fitting parameter. 

In summary, we have shown that it is possible to relate 
the second moment of the g factor of metallic nanopar- 
ticles with strong spin-orbit coupling to other spectral 
statistics of energy levels without resorting to any micro- 
scopic parameter. Our results also show that it is possi- 
ble to estimate the spin and orbital contributions to the 
fluctuations of the g factor by comparing data taken from 
nanoparticles doped and undoped with a heavy-element 
metal. We suggest that a fitting of the data to a univer- 
sal, dimensionless level velocity correlation function may 
provide an additional test of the applicability of random 
matrix theory to these systems and allow us to extract 
information about the intragrain disorder. 

This work was supported in part by NSF grants No. 
DMR 02-37296, No. DMR 04-39026 (L.I.G.) and No. 
CCF 0523603 (E.R.M.). C.H.L. acknowledges partial 
support in Brazil from CNPq, Instituto do Milenio de 
Nanociencias, and FAPERJ. E.R.M. acknowledges par- 
tial support from the Interdisciplinary Information Sci- 
ence and Technology Laboratory (I 2 Lab) at UCF. We 
are grateful to Y. Fyodorov, J. Petta, and F. von Oppen 
for enlightening discussions. L.I.G. and E.R.M. thank 
the Instituto de Fisica at UERJ, Brazil, and the Aspen 
Center for Physics for the hospitality. 



1 D. G. Salinas, S. Gueron, D. C. Ralph, C. T. Black, and 
M. Tinkham, Phys. Rev. B 60, 6137 (1990). 

2 D. Davidovic and M. Tinkham, Phys. Rev. Lett. 83, 1644 
(1999). 

3 P. W. Brouwer, X. Waintal, and B. I. Halperin, Phys. Rev. 
Lett. 85, 369 (2000). 

4 K. A. Matveev, L. I. Glazman, and A. I. Larkin, Phys. 
Rev. Lett. 85, 2789 (2000). 

5 R. A. Serota, Solid State. Comm. 117, 605 (2001). 

6 J. R. Petta and D. C. Ralph, Phys. Rev. Lett. 87, 266801 

(2001) . 

7 J. R. Petta and D. C. Ralph, Phys. Rev. Lett. 89, 156802 

(2002) . 

8 A. Szafer and B. L. Altshuler, Phys. Rev. Lett. 70, 587 
(1993). 

9 An analogous relation was derived for the variance of the 
persistent current of a mesoscopic ring in the presence of 
spin-orbit coupling by V. E. Kravtsov and M. R. Zirnbauer, 
Phys. Rev. B 46, 4332 (1992). 

The magnetic moment operator is antisymmetric with re- 
spect to time reversal. Therefore, eigenvectors related by 



time reversal have eigenvalues equal in amplitude but with 
opposite signs. 

C. E. Porter, in Statistical Theories of Spectra: Fluctua- 
tions, edited by C. E. Porter (Academic Press, New York, 
1965). 

S. Adam, M. L. Polianski, X. Waintal, and P. W. Brouwer, 
Phys. Rev. B 66, 195412 (2002). 

M. L. Mehta, Random Matrices, 3rd edition (Academic 
Press, Amsterdam, 2004). 

At B = 0, due to Kramers degeneracy, 5 is equal to half of 
the mean level spacing. 

Y. V. Fyodorov and H.-J. Sommers, Z. Phys. B 99, 123 
(1995). 

F. von Oppen, Phys. Rev. E 51, 2647, (1995). 

B. D. Simons and B. L. Altshuler, Phys. Rev. Lett. 70, 

4063 (1993). 

B. D. Simons and B. L. Altshuler, Phys. Rev. B 48, 5422 
(1993). 

I. E. Smolyarenko and B. D. Simons, J. Phys. A 36, 3551 
(2003). 



